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1. Introduction 

These lectures arc an introduction to the subject of sum rules in hadronic 
interactions as understood at present within the framework of quantum 
chromodynamics (QCD). The literature on the subject is extensive and 
various applications are quoted in other lectures at this school. I have 
chosen to concentrate on the foundations of the subject rather than giving 
an exhaustive list of results obtained from sum rules. The aim is to give 
you a critical overview of the subject so that you can judge by yourselves 
on how seriously to take the results in the literature, and also to foresee 
which possible new applications could be made. One often encounters two 
extreme attitudes among theorists on results obtained using QCD sum 
rules. One is the optimistic belief that the result can be trusted at the few 
percent level, when in fact there are a host of assumptions in the derivation; 
the other is the pessimistic attitude which pretends that all QCD sum rules 
are "wrong" because some calculation somebody made somewhere gives an 
incorrect result. I want to show you that the existence of sum rules in many 
cases follows rigorously from QCD. Often the difficulty is not in the sum 
rule itself, but rather on how to use it in practice, when one is dealing with 
a channel where there is very little direct experimental information. We 
shall sec that it is important to check that the matching of the input low 
energy hadronic ansatz to the continuum described by perturbative QCD 
(pQCD) satisfies global duality tests. 

The lectures are organized as follows. I shall start in section 2 with an 
introduction to sum rules which were derived prior to the development of 
QCD. This section is mostly descriptive and assumes some familiarity with 
topics like chiral perturbation theory which arc covered in other lectures [1, 
2] at this school, and with dispersion relations. The reader who has never 
heard of dispersion relations should first read the beginning of section 3 and 
section 3.1, and then go back to section 2. Section 4 contains quite a few 
technical details. It should be useful to those who want to learn how to use 
QCD sum rules in practice. Section 5 gives a qualitative introduction to 
non-perturbative power corrections. It should be helpful as an introduction 
to more technical papers quoted there. The two examples discussed in 
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section 6 arc rather topical subjects, much under discussion at the moment, 
and illustrate how to use the material covered in the previous sections. 

There are a lot of topics on sum rules which are not mentioned in these 
introductory lectures like for example, applications to baryons, to heavy 
quark systems and gluonium. The curious reader should complement these 
lectures with extra reading. There are two references which contain a 
lot of material. One is the book by Stephan Narison [3]; the other is 
M.A. Shifman's selection of some of the original articles with his comments 
in ref. [4]. 



2. Sum Rules prior to QCD 



Sum rules in hadron physics have a long history. We shall be concerned 
in this section with various sum rules which were suggested before QCD 
was established as the theory of the strong interactions. I want to do this 
for two reasons: first, because they are interesting sum rules per se and 
second, because they illustrate very clearly the basic ingredients which go 
into the derivation of sum rules in general. The DrcU Hearn sum rule and 
the Adler-Weisberger sum rule are expected to be true in QCD. As we 
shall see later in section 5.3, the Weinberg sum rules are in fact theorems 
of QCD in the chiral limit where the light quark masses are neglected. 



2.1. The Drell-Hearn Sum Rule 



The sum rule in question [5] relates the proton anomalous magnetic moment 
^2(0) to an integral over the energy of the difference of total photoabsorp- 
tion polarized cross-sections 



[^^(0)]' =^r^ Ki>(-) - ■ (2-1) 



■th 



Here, ctTT(i/) [cr^|,(i/)] denote the total cross-section for the absorption of 

a circularly polarized photon of laboratory energy by a proton polarized 
with its spin parallel (antiparallel) to the photon spin, Mp is the proton 
mass and a ~ 1/137 the fine structure constant. The integral starts at the 
inelastic threshold i/^j^ = TO^/2Mp. 

The three ingredients which are needed to derive this sum rule are: 
dispersion relations for forward Compton-scattering, the optical theorem 
which relates total cross-sections to the imaginary part of forward scatter- 
ing amplitudes, and a low-energy theorem which fixes the value of the real 
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forward scattering amplitude at zero energy. We shall see these three ingre- 
dients appearing systematically in the derivation of all the sum rules which 
we shall consider. However, in many of the QCD sum rules which we shall 
later discuss the low energy theorem, will be replaced by a high energy the- 
orem instead, or rather a short-distance behaviour property which exploits 
the fact that two-point functions in QCD are known perturbatively at high 
euclidean momentum transfer due to the asymptotic freedom property of 
QCD. 

The forward Compton-scattering amplitude illustrated by the diagram 
in Fig. 1, 




P P 

Fig. 1. Forwcird Compton scattering of photons 7(g) on protons p. 

where p ■ q = Mpi/, depends on two scalar invariant amplitudes of the 
squared energy v'^: 

= h{v^)e' ■ e + uf2{u^)ia ■ e' x e. (2.2) 

The spin-flip amplitude obeys a dispersion relation which is expected to 
have no subtractions i.e., 

Re/2(z/2) = PP / -Im/2(z>2). (2.3) 

[Section 3 will be dedicated to the study of dispersion relations in general. 
They are the basics of sum rule derivations.] The optical theorem relates 
Im/2(i'^) to total cross-sections: 

ST,^h{v-')=cjl\,{u)-a\l{u), (2.4) 

and the low energy-theorem relates Re/2(0) to the proton anomalous mag- 
netic moment squared: 

/2(0) = -^^[F2(0)f . (2.5) 
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The optical theorem is a general property of unitarity in quantum field 
theory. The low-energy theorem follows from chiral symmetry properties 
of the underlying QCD Lagrangian [See the lectures of Pich [1] and Manohar [2] 
in this school.] The dispersion relation for forward scattering amplitudes is 
also a general property of quantum field theory [6]. What is still lacking 
to promote this sum rule to the status of a QCD theorem is the proof 
that Im/2(j^^) — > when v ^ oo. At present this accepted behaviour is 
only supported by theoretical strong interaction wisdom, like Regge phe- 
nomenology. 

The Drell-Hearn sum rule has been checked in perturbation theory in 
QED-like theories to the first non-trivial order. One may ask if this sum 
rule could be used for practical calculations of the anomalous magnetic 
moment of the electron. The answer is that it is not a competitive method, 
simply because F2(0) in the l.h.s. appears squared which obliges one to 
calculate the integral in the r.h.s. to very high orders in powers of the fine 
structure constant. 

The Drell-Hearn sum rule is also useful in the phenomenological study 
of deep inelastic scattering of polarized leptons on polarized protons [See 
Manohar's lectures [2] in this school.], where it provides a constraint on the limit 
of real photoabsorption. It plays an important role as well in the calculation 
of the so called polarizability contribution to the hyperfine splitting in the 
hydrogen atom [See e.g. ref. [7] and references therein.] 

2.2. The Adler-Weisberger Sum Rule 

The sum rule in question [8,9] relates the axial coupling constant of 
/3-decay to integrals of pion-nucleon total cross-sections: 



1 a^-^Jl 



'tot 



'^W-^t'^t^M , (2.6) 



where ~ 92MeV denotes the pion decay coupling constant. The refer- 
ence amplitude in this case is the forward pion-nucleon scattering ampli- 
tude illustrated in Fig. 2 
which has two invariant amplitudes: 

T,j{u) = 5,jT^+Hu) + \ [r,,T,-] ^T(-)(z.), (2.7) 

where z,j are isospin indices. The basic dispersion relation is the one 
obeyed by the isospin-odd amplitude T'^~\v), and as in the case of the 

Compton scattering spin-flip amplitude, it is expected that T'^~^{v) obeys 
an unsubtracted dispersion relation. The optical theorem relates ImT'^~) {v) 
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_7r(q) 



7U(q) 




N(p) N(p) 

Fig. 2. Forward pion-nucleon scattering amplitude. 



to the difference of 7r~-p and 7r+-p total cross-sections. On the other hand, 
the calculation of the low-cncrgy behaviour of the amplitude ReT'^~)(j^) 
can be done using baryon chiral perturbation theory [2,1]. The relevant 
Feynman diagrams are those of the Born approximation shown in Fig. 3 



7t(q) 7t(q) 7U(q) 




N(p) N(p) 



Fig. 3. Forward pion-nucleon scattering at the Born approximation. 



with the result 

limReT(-)(z/) = l-g\. 



(2.8) 



The seagull graph contributes the factor 1, each one of the other graphs a 
factor —\g\- 

Sometimes, the Adler-Weisberger sum rule is also written in the follow- 
ing way 



1_ J_ = i^ 

5a T^gln^ 



Jo V 



(2.9) 



Both forms are equivalent, because of the Goldberger-Treiman relation [10] 

SIaMn = ffTrNN/TT- (2.10) 

The Goldberger-Treiman relation is also a low-energy theorem of QCD. 

I would like to come back to the question of an unsubtractcd dispersion 
relation for the isospin- odd amplitude T'^-\v). This, of course depends 
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on the high-energy behaviour of the underlying theory. One can consider 
as a toy model the linear sigma model [11] for constituent quarks. It has 
been shown [12,13] that in perturbation theory, this model does not require 
subtractions in the corresponding Adler-Weisberger sum rule and yields, 
at the one-loop level, the result 



51 = 1 - log ttI + 0[l) , (2.11) 




where Mq denotes the mass of the constituent quark in this model. On 
the other hand, if one considers an effective constituent chiral quark model 
like the Georgi-Manohar model [14], then it can be shown [15] that the 
corresponding isospin-odd amplitude T*-^^ does not obey an unsubtracted 
dispersion relation. Beyond the tree level, ReT(~)(j/) cannot be obtained 
from a dispersive integral; and indeed, the calculation at the one loop level 
reproduces the result in eq. (2.11) but with a cut-off replacing the sigma 
mass. It is simply that in this case the effective low-energy theory does 
not have the high-energy behaviour of the underlying theory. 

2.3. The Weinberg Sum Rules 

There is a correlation function which is particularly sensitive to properties 
of chiral symmetry breaking, namely the two-point function 

^ZiiQ) = '^ij rf*a;e'«-(0|T {L<^{x)R''{Q)^) |0), (2.12) 

with currents 

= J(x)7^i(l-75)w(a;) and R^" = d{x)Y\{l + l5)u{x) . (2.13) 

In the chiral limit where the light quark masses are set to zero, this two- 
point function only depends on one invariant amplitude (Q^ = —q^ > for 
spacelike) 

Kr{Q^) = W - 9''''q^)nLR{Q^). (2.14) 

The function Hlr{Q'^) vanishes order by order in perturbation theory and 
is an order parameter of the spontaneous breakdown of chiral symmetry 
(SxSB) for all values of the momentum transfer. As illustrated in Fig. 4 
the two point function in eq. (2.12) is the difference of a vector current cor- 
relation function and an axial-vector current correlation function. Chiral 
symmetry identifies these two correlation functions. 

Two-point functions in general obey dispersion relations. [Scc the discus- 
sion in section 3.] In 1967, and therefore prior to the development of QCD, 
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R = VA i<v 



Fig. 4. The LR two-point function is the difference of VV and AA two-point functions. 

Weinberg made an educated guess about the convergence properties of the 

LR two point function, with the suggestion that the following two sum 
rules may hold in the underlying theory of the strong interactions [16]: 

/■oo /'OO 

/ dt\mRLR{t) = dt [ImUvit) - Imn^Ct)] = 0, (2.15) 
Jo Jo 

and 

/•oo />oo 

/ dttlmn.LR{t)= j rftt[Imny(t) -lmnA(t)] = 0. (2.16) 
Jo Jo 

They are commonly referred to as the 1st Weinberg sum rule and the 2nd 

Weinberg sum rule, and they arc examples of the so called superconver- 
gent sum rules for the following reason. Let us assume that the function 
^lr{Q^) obeys an unsubtracted dispersion relation: 
dt 

TIlr{Q^) = ^^ImULRit). (2.17) 

The existence of the two sum rules above implies a stronger convergence 
behaviour for the function IIlr{Q^) in the deep euclidean region than 
the minimum required for the unsubtracted dispersion relation to hold. 
Expanding j:fQ2 — — -^t-^ + ■ ■ ■ , we observe that the 1st Weinberg 
sum rule corresponds to the constraint 

lim Q^nLR{Q^)^0, (2.18) 

and the 2nd Weinberg sum rule to the further constraint 

lim g^nifl(Q2)^o. (2.19) 

Q^^oo 

It is quite remarkable that, as we shall see in section 5.3, both constraints 
are now QCD properties; and in fact, the 1st Weinberg sum rule still holds 
in the presence of finite light quark masses [17]. 

I would like to make some comments on these sum rules. 

i) Weinberg in his original paper [16] assumed that only the lowest nar- 
row resonant state in the vector and axial-vector spectral functions: 

-ImUv {t) = f^M^6{t - M2 ) + . . . (2.20) 
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and 

-ImUAit) = fj{t) + flMl5{t - Ml) + • • • , (2.21) 

TT 

contribute significantly to the sum rules. The first term in the axial-vector 
spectral function in eq. (2.21) is the contribution from the massless pion. 
Thus the 1st and 2nd Weinberg sum rules in eqs.(2.15) and (2.16) constrain 
the couplings and masses of the narrow resonances as follows 

fyMt.-flMX = 0. 

It is interesting that the original Weinberg's ansatz in eqs. (2.20) and (2.21) 
can now be justified within the framework of the large-A''c limit of QCD. 
Indeed, in this limit, the vector and axial-vector spectral functions are 
each a sum of an infinite number of narrow sates. The dots in the spectral 
functions in eqs. (2.20) and (2.21) implicitly assume that the sum of the rest 
of the narrow states is already globally dual to the onset of the perturbative 
continuum, 

hmUvit) = f^M^6{t - M2) + ^le{t - so)[l + • • •], (2.23) 



and 



hmUAit) = fj{t) + f\Ml5{t - Ml) + ^\0{t - so)[l + •••], 

(2.24) 

where the dots indicate now gluonic corrections. The vector and axial- 
vector spectral functions of the pQCD continuum are the same in the chiral 
limit. 

ii) In order to make a "prediction" for the axial mass Ma, Weinberg in 
his original paper also makes the assumption that 



/^ = 2^, (2.25) 



which was called at the time the KSFR relation [18,19]. This, together 
with eqs. (2.22) leads to the prediction Ma = sPlMy which is satisfied by 
the central values of the presently known masses within an acciiracy of less 
than ~ 12%. It can be shown [20] that the KSFR-relation follows from the 
assumption that both the electromagnetic pion form factor and the axial 
form factor in tt ^ c-Vel obey unsubtracted dispersion relations, plus the 
dynamical input of resonance dominance. 
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iii) Inverse moments of the difference of vector and axial— vector spectral 
functions with the pion pole removed, i.e. integrals like 



withp = 1, 2, 3, • • • , correspond to non-local order parameters which govern 
the couplings of local operators of higher and higher dimensions in the low 
energy effective chiral Lagrangian [See Pich's lectures [1] in this school.]. For 
example, the first inverse moment is related to the 0{p^) coupling constant 
I/io as follows [21] 



where the tilde in Imll indicates that the pion pole has been removed. 
Using the ansatz in eqs. (2.23) and (2.24) and the two Weinberg sum rules 
in eq. (2.22) one finds 



which provides a good estimate [20] of the iio coupling constant. 

iv) One of the parameters which characterize possible deviations from 
the Standard Model predictions in the electroweak sector is the so called S- 
parameter [See the lectures of Treille [22] and Chivukula [23] in this school.] . In the 
unitary gauge, S measures the strength of an anomalous Wll,/ B^'^ coupling. 
This is the SU (2) l x SU (2)i^ coupling analogous to the term proportional to 
Lio in QCD. It has been argued that if the underlying theory of electroweak 
breaking is a vector-like gauge theory of the QCD type, the sign of the 
Lio-like coupling must be negative, precisely as observed in QCD. In the 
electroweak sector, this fact is infirmed by experimental observation and 
constitutes at present a serious phcnomcnological obstacle to technicolour- 
like formulations of electroweak symmetry breaking. The relation of the 
Lio coupling to the relative size of local order parameters versus in the 
large-TVc limit of QCD-like theories has been recently discussed in ref. [24] . 

2.4- The Electromagnetic Pion Mass Difference 

The most remarkable application of the Weinberg sum rules is the calcula- 
tion of the TT"*" — 7r° electromagnetic mass difference [25] . In the chiral limit 
and to lowest order in the electromagnetic interactions, there appears a new 
term with no derivatives in the low-energy effective chiral Lagrangian [26] 




(2.26) 




(2.27) 




(2.28) 



jC^f[ = --- + e^Ctr{QRUQLU^) 



(2.29) 
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where U is the matrix field which collects the octet of pseudoscalar 

Goldstone fields [See Pich's lectures [1] in this school.] and Qr = Ql = 

diag[2/3, — 1/3, — 1/3] the right- and left-charges associated with the elec- 
tromagnetic couplings of the light quarks. Expanding U in powers of the 
pseudoscalar fields there appear quadratic terms like 

£eff = • • • - 2e2C7-l(7r+7r- + K+R-) + ■■■ , (2.30) 

/ TT 

showing that, in the presence of electromagnetic interactions, the charged 
pion and kaon fields become massive: 

2e^C 

(m^+ - m^o)|EM = (m^+ - m^o)|EM = -J^- (2-31) 

In fact, the main contribution to the physical tt^-tt" mass difference is of 
electromagnetic origin because quark masses do not contribute significantly 
to the 7r+ tt" mass splitting. 

It can be formally shown [See e.g. refs. [27,28] and references therein.] that 
the coupling constant C of the effective term in eq. (2.29) which results 
from the integration of a virtual photon of cuclidcan momentum squared 
in the presence of the strong interactions in the chiral limit, is given by 
an integral of the correlation function IlLRiQ"^) defined in eqs. (2.12) and 
(2.14): 



87r2 4 



/ dQ^ Q^IIlr{Q^). (2.32) 

JQ 



The constraints in eqs. (2.18) and (2.19) guarantee the convergence of the 
integral in the ultraviolet. Inserting the resonance dominance ansatz in 
eqs. (2.23) and (2.24) for the vector and axial-vector spectral functions in 
the dispersion integral representation in eq. (2.17) results in 



^LR{Ci ) - - ^^T^ + ^ + (2-33) 



With this expression inserted in the integrand of the r.h.s. in eq. (2.32) and 
using the 1st and 2nd Weinberg sum rules in eqs. (2.22) one easily obtains 

(rn^.-m^c)|KM = -^^^^^log^. (2.34) 



The authors of ref. [25] also use Ma = \f2Mv which results in 
a3 

TT 2 



(m2+ - m^o)|EM = r o log 2- (2.35) 
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Numerically, this corresponds to a mass difference AmirlEM — 5.2 MeV to 
be compared with the experimental mass difference —fn^o = 4.6 MeV. 

Equation (2.32) can nowadays be viewed as a rigorous QCD sum rule 
in the chiral limit. The sum rule relates the low energy constant (to^+ — 
rn^o ) I EM to an integral of the real part of a two-point function. The integral 
runs over all possible values of the euclidean momenta and therefore it 
requires, a priori, the knowledge of the two point function at all distances. 
I want to stress the fact that it is not enough to know the very low 
behaviour of the \iLR{Q'^) correlation function to obtain an approximation 
to the r.h.s. integral in eq (2.32). The low behaviour of this function is 
governed by chiral perturbation theory and the first two terms are known 

nLR{Q'') = -^-ALw + 0{Q^); (2.36) 

but this behaviour cannot be extrapolated to the ultraviolet region of the 
integral since it would lead to a divergent result. By contrast, as we shall 
further discuss in section 5.3, the full resonance dominance ansatz in 
eq. (2.33) with the Weinberg sum rule constraints in eqs. (2.22) has both 
the long-distance behaviour and the short-distance behaviour of QCD. 



3. Dispersion Relations 

In these lectures we shall be all the time considering two-point functions 
of local operators i.e., Fourier transforms of the vacuum expectation value 
of the time ordered product of a local current J{x) times its hermitian 
conjugate 

n(g2) = i j rfW-^(0|T { J(a;) J(0)t} |0). (3.1) 

In most applications the current J{x) is one of the Noether currents asso- 
ciated with global gauge transformations of flavour degrees of freedom, like 
a vector current q{x)"/'^q{x), or an axial vector current q{x)"/^"f5q{x): but 
one could also consider gauge invariant local operators of gluon fields like 
tTG^'^{x)Gfj,v{x), or composite operators like q{x)Tq{x)q{x)rq{x). It has 
been shown by Kallen and Lehmann [29,30] quite a long time ago that two- 
point functions obey dispersion relations. The dispersion relation follows 
from the analyticity properties of 11 (g'^) as a complex function of q"^, the 
only energy-momentum invariant which appears in a two-point function. 
In full generality n(g^) is an analytic function in the complex g^-plane but 
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for a cut in the real axis < < oo. It then follows that 

n((7^) = / dt- 5 -lmn(i) + a + b?^ + . . . , (3.2) 

where the degree of the arbitrary polynomial in the r.h.s. depends on 
the convergence properties of lmn{t) when t — > oo. The interest of this 
representation is that -!-Imn(t) in the integrand, which is usually called 
the spectral function, is a physical cross-section. With J{x) a current with 
specific quantum numbers, the spectral function is then directly related 
to the total cross section for the production from the vacuum of hadronic 
states with those quantum numbers. 

For example, with J the electromagnetic hadronic current of light 
quarks, 

J^'{x) = ^uix)^^'u{x) - ^dixj^^dix) - \s{x)'y^'s{x) (3.3) 

o o o 

the relation to the total e+e" annihilation cross-section into hadrons is 
<^{q )e+e-^hadrons = ^ -Imn^Mlg ), (3.4) 



and 



-3^(g)g2 ilmnEM(g') = 

TT 



r 



(3.5) 



where the sum is extended to all possible physical states, on-shell states, 
with an integration over their corresponding phase space understood. A 
pictorial representation of these equations is shown in Fig. 5 below. 




Fig. 5. Representation of the hadronic spectral function in eq. (3.5). 



The lowest possible state in this case is a two-pion state. The function 

Hem (9^) is analytic in the complex -plane but for a cut in the real axis: 
4to^ < q"^ < 00, as illustrated in Fig. 6. 
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Physical Region 



Complex q 2 -Plane 



Fig. 6. Representation of the complex q'^-plane. 



I want to show first how these properties appear explicitly in a simple 
case which you should be familiar with, and then we shall discuss a general 
proof. 

3.1. Vacuum Polarization due to Massive Fermions 

Let us consider QED with a massive fcrmion. We want to calculate the 
lowest order vacuum polarization contribution from this fermion, i.e. the 
one-loop Feynman diagram in Fig. 7. 



Combining the two fermion propagators with a Feynman parameter x, 



and after integration over the four-momenta k in the loop [l am assuming 
that a gauge invariant regularization of this integral, like e.g. dimensional regulariza- 
tion, has been made.], the one-particlc irreducible part of this diagram, when 
subtracted at = if we choose on-shell renormalization, leads to the 




Fig. 7. One fermion loop contribution to vacuum polarization. 




2 ' 



(3.6) 



18 



E. de Rafael 



renormalized self-energy function: 

HrC?^) = - ^ ^ dx2x{l-x)\og(^- .^ x{l - x)j , (3.7) 

With the change of variables y = 1 — 2x, and using the fact that the 
resulting integral is symmetric when y — > —y, we get 

U^{q^) = - [\y{l-y^)\og' 

Jo 

tegrating by parts this equ; 
results in a rational integral: 



Integrating by parts this equation using the identity: I — y^ = ■§-{y— ky^) 



(3.8) 



n^(-') = ^l^«M,-lv'} ,^,_J_,..^_,^ , (3.9) 

and with a new change of variables: t = wc finally get a rather 

interesting representation of the renormalized self energy 



nH((?') 



a f°° dt 1 1 / 2rrf\ I 4m^ 
^ JAm^ t-q^ -ie ~T~ J V ~r' ^ ■ ' 

The reason why this representation is interesting is that it is in fact a 
dispersion relation! We have succeeded in rewriting the initial Feynman 
parametric representation in eq. (3.7) as a dispersion relation by simple 
changes of variables. Using the identity 

\ -=PP-^+i7r(5(i-g2), (3.11) 

we immediately see that 

llmn(t) = + ^) f^Oit - Am-). (3.12) 

Equation (3.10) is a particular case of the general dispersion relation writ- 
ten in eq. (3.2), when the arbitrary polynomial is just a constant, and we 
have got rid of the constant because the on-shell renormalized IIr is defined 
as 

n^(92) = n(92) - n(o) = / ^ \ . -imn(t). (3.13) 

^0 t t — q'^ — te TT 
It is perhaps worth insisting on the fact that asymptotically 
1 a 1 

lim -lmn(t) ^ ; (3.14) 

i.e. the electromagnetic spectral function goes to a constant. In fact, it is 
this constant which fixes the value of the lowest order contribution to the 
/3-function associated with charge renormalization in QED. 
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3.2. General Proof of Dispersion Relations 

We shall now sketch a proof of the dispersion relation property for two- 
point functions in full generality. The key of the proof lies in the definition 
of the time ordered product implicit in eq. (3.1): 

T(J(x) J(0)^ = e{x)J{x)J{Q)^ + 0{-x)J{O)^J{x), (3.15) 

and the use of translation invariance. The function 9{x) is the Heaviside 
function: 9{x) = 1 if xo > and ^(a;) = if a;o < 0, which has the integral 
representation 

First we insert a complete set of states |r)(r| between the two cur- 
rents in the T-product definition. This leads to matrix elements of the 
type (0| J(x)|r) to which we apply translation invariance: 

(0|J(.T)|r) = {{)\U-^UJ{t)U-^U\T), (3.17) 

where U is the unitary operator induced by translations in space-time: 

U{a)J{x)U-^{a) = J{x + a) and Z^(a)|r) = e'f^"|r), (3.18) 

and pt denotes the sum of the energy-momenta of all the particles which 

define the state |r). The choice a = —x factors out the a;-dependence of 
the matrix element into an exponential: 

(0| J(x)|r) e-'P^-=^ (0| J(0)|r). (3.19) 

All the particles in the state |r) are on-shell. This constrains the total 
energy-momentum pr to be a time-like vector: p^ = t with t > 0. With 
these constrains on pr we can insert the identity 

dt [ cf^p 9{p) Sip"^ - t) {p - Pr) = 1 (3.20) 







inside the sum J2r o^^'" complete set of states. Interchanging the order 

of sum over F and integration over t and p, there appears naturally the 
definition of the spectral function associated with the J-operator 

^(0|J(0)|F)(F|J(0)t|0) (27r)45W(p-pr) ^ 27rp(/). (3.21) 
r 

The spectral function p{p^) is a scalar function of the Lorentz invariant p"^ 
and the masses of the particles in the states |F) only. By construction it is 
a real function and non-negative 

p{p'r = p(p2) > 0. (3.22) 
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We can now rewrite the two-point function in eq. (3.1) as follows 
n{q^) = J d'^x e'«-^ dtp{t) X 



/ 



^ [ie{x) e-^f-e(p)^(/ -t)+ ie{-x) e'P-^e{p)6ip^ - t)] . 

(3.23) 

Here, you may recognize the familiar functions of free field theory 

^^^"^^ = / (S^ e-^^'-^(p)5(p2 - t) (3.24) 
and 

^~{x)= I 4^e^'' ''Kp)5{p'-t) 



d^p 
(27r)3 

'^'P .e-^P■^e{-p)5{p^-t)■ 



{2,.2h) 



(27r)3 

and therefore the Feynman propagator function 
Ap(a:;t) = ie{x)A'^ {x\t) + ie{-x)A- {x;t) 

d'^p e-*P-^ (3-26) 



/ 



[iTvyt-p^-ie 



where the last expression can be obtained using the representation in 
eq. (3.16) of the 6 function [Scc e.g. rcf. [,3l] pp. 50-51.]. The two-point 
function n(g^) appears then to be the Fourier transform of a scalar free 
field propagating with an arbitrary mass squared t weighted by the spectral 
function density p{t) and integrated over all possible values of t, 

n(g2) = y d^j^gig-x^ dtp{t)A-p{x;t). (3.27) 

Integrating over x and p results finally in the wanted representation 

1 

n(g2) = / dtp{t)- ^— . (3.28) 

Jo t-q^ -le 

With 

Il{q^) = Ren(g^) + ilmU{q^), (3.29) 
and the use of the identity in eq. (3.11), it follows that 

p{t) = -lmn(t), (3.30) 
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which identifies the spectral function with the imaginary part of the two- 
point function. 

Notice that the formal manipulations above avoid the question of con- 
vergence of the principal value integral 



The convergence of the integral in the ultraviolet limit {t ^ oo) depends on 
the behaviour of the spectral function at large t-values. When doing above 
the exchange of sum over F and integrations we have implicitly assiimcd 
good convergence properties; but in general the product of the distributions 
6{x) and dt p{t) /S.'^ {x;t) may not be a well-defined distribution. The 
ambiguity manifests by the presence of an arbitrary polynomial in in 
the r.h.s. of the PP-integral 



Notice that the coefficients of the arbitrary polynomial have no disconti- 
nuities; in other words, the ambiguity of short distance behaviour refiects 
only in the evaluation of the real part of the two-point function, not in 
the imaginary part. The physical meaning of these coefficients depends of 
course on the choice of the local operator J{x) in the two-point function. 
In some cases the coefiacients in question are fixed by low-energy theorems; 
e.g. if n(0) is known, we can trade the constant a in eq. (3.32) for 11(0): 



In other cases the constants are absorbed by renormalization constants. In 
general, it is always possible to get rid of the polynomial terms by taking 
an appropriate number of derivatives with respect to q^. Various examples 
will be discussed later. 

3.3. QCD Moment Sum Rules 

In QCD the number of derivatives which is required to obtain a well defined 
two-point function is fixed by the asymptotic freedom property of the the- 
ory. For a gauge invariant local operator J{x), the asymptotic behaviour 
of the associated two-point function is of the type 




(3.31) 




(3.32) 




(3.33) 




(3.34) 
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with A and ai calculable cocfScicnts, and d a known integer d = 0, 1, 2, . . . , 
depending on the dimensions of the operator J{x). It is then sufScient 
to take d + 1 derivatives with respect to <f' to get rid of the arbitrary 
polynomial and obtain a convergent integral. The functions defined by the 
moment integrals (Q^ = —q^): 



POO 

^ Jo (t + W+i \t + Qy n 

(3.35) 

for m > (i + 1 arc then well defined functions calculable in pQCD at sufla- 
ciently large Q^-values. 

3.4- The Adler Function 

An example of a QCD moment sum rule is the Adler function correspond- 
ing to the two-point function associated with the vector-isovector quark 
current 

^1=1 {x) = ^ [: u{x)')i'u{x) : - : d{x)'yi'd{x) :] . (3.36) 
Since the current is conserved, the two-point function 

nr:i(a) = ^ I d''xe^'> ^Q\T{Vt,{x)V{^M} |0> ^ ^ 

J (3.37) 

has one invariant function ni=i(<7^) only. Power counting tells us that in 
this case the asymptotic behaviour of the corresponding spectral function 
goes as a constant i.e. the power d in the r.h.s. of cq. (3.34) is d = 0. 
Only one derivative is required to have a well defined correlation function. 
The Adler function is defined as the logarithmic derivative of n(g^)i=i 
(Q^ = -q^, with > for g^-spacelike), 

AQ') - -Q'^^^P- (3-38) 

This function, for suSiciently large values of is calculable in pQCD as 

a power series in the running coupling constant. On the other hand the 
vector-isovector spectral function ilmni=i(i) is accessible to experimental 
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determination from e+e annihilation into hadrons and from the r hadronic 
decays. We then have the sum rule 



When trying to confront this sum rule with experiment, there appears 
the problem that the integrand in the r.h.s. is only known experimentally 
from threshold up to finite values of t. This brings in a question of matching 
whatever is known about the low energy hadronic spectral function with 
its asymptotic behaviour as predicted by pQCD. We shall later on come 
back to this important question, after the discussion of non-perturbative 
power corrections. 

The Adlcr function has been calculated in perturbation theory to the 
four-loop level [32-36]. In the MS'-scheme and with = a-j^^Q"^), the 
numerical evaluation of the exact perturbative calculations for Nc = 3 and 
for three light flavours n/ = 3 gives the result 



4. Types of Two Point Function Sum Rules 

There are several types of sum rules which have been discussed in the 
literature. They differ on the way that the analyticity properties of two- 
point functions are exploited and/or on various limits that one can define 
on moment sum rules. We explain the most common types of sum rules 
below. 

4-1- Laplace Transform Sum Rules 

The Laplace transform of a spectral function: 



is a calculable function in pQCD provided that the variable a, which has 
the dimensions of an inverse mass squared, is sufficiently small. The inter- 
esting point about this type of sum rules is the presence of the exponential 
factor in the integrand which gives a predominant weight to the low-energy 
component of the hadronic spectral function compared to the power weight 
of a moment sum rule. Sum rules of this type were first proposed by the 
ITEP group [37] who called them Borel sum rules. The Laplace transform 




(3.39) 




(3.40) 




(4.1) 



24 



E. de Rafael 



in (4.1) results from a precise limit of the moment sum rules discussed in 
section 3.3; it is the limit where the number of derivatives of the invari- 
ant two-point function n(g^) with respect to Q'^ and Q"^ = —(f' both go to 

infinity with their ratio ^ = = i fixed. More precisely, the limit is 
obtained by applying the operator 

(iv-i)!^^^ (pq^ 



L = lim lim I w 



(4.2) 



to the first well-defined moment of n(g^), 
behaviour defined in eq. (3.34)]: 



d is the asymptotic power 



1l 

(7 



(-1) 



Qd+1 



L 



°° N{N + 1) 

(t+WTT 



n(g^) 

(N + d+l) 



(4.3) 



1 + 



N 



-lmn(i). 



In the limit which we are considering, the first fraction in the integrand 
goes to cr'^+^ and the second fraction to e~*^, with the result that 



1l 

a 



(-1) 



d+1 



Qd+1 



(5Q2) 



2^d+l 



n(g2) 



d+1 



dte-*'^ -linU{t). (4.4) 



In practice, the question which then appears is how to calculate the 
function A4 (cr) without having to take the complicated limits which define 
the operator L in eq. (4.2). This is possible because L is nothing but an al- 
gebraic form of the Laplace inversion operator. Once this is recognized [38] , 
it is easy to go directly from the moment sum rule 



ad+l i-oo 



{d + iy . 1 

(t + Q2) 



2\d+2 



lmU{t) 



(4.5) 



to the Laplace transform in (4.1). Indeed, the power weight in the integrand 
of eq. (4.5) can be written as a Laplace transform with respect to the 
variable 



(d+1)! 



{t + Q2)d+2 - 

and the action of L brings back the inverse Laplace transform i.e., 
((d + 1)! 



(i+Q2)d+2 



(4.6) 



(4.7) 
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This provides the key to calculate the Laplace transform M{a) once we 
know n(g^). The results of analytic pQCD calculations of self-energy func- 
tions like n((j'^) are combinations of terms which are typically of the form: 



V^, ...}. (4.8) 

The trick to go from Ti{q^) to M{(t) is to write each of these terms as a 

Laplace transform with respect to tli(^ variable . Since the operator L 
acts in fact as an inverse Laplace transform one then finds e.g. that 



(Q 



1 , 

2\a+l „2 



r(a + 1) 



-(-logc7iv2+^(a + l)), (4.9) 



where ^(z) = ^ logr(2). [Other useful formulae like this one can be found e.g. in 
ref. [38].] 

There arc a number of interesting properties of the Laplace transform 
which are useful to know when doing phenomenological applications and 
therefore I want to discuss them next. In sum rule applications one often 
considers the ratio [39] 

At very small values of a, the function 1Z{a) has the free-field behaviour 
7^(a)^^(l + ...), (4.11) 

(7 

where the dots indicate perturbative gluonic corrections. Recall that d is 
the asymptotic power behaviour in t of the corresponding spectral function 
as shown in eq. (3.34). (For example d = in the case where J{x) is the 
vector— isovector current.) This dependence in d is the only reminiscence left 
from the subtractions needed in the dispersion relation for the initial two- 
point function. At very large a the behaviour of 7?.((j) will be dominated 
by its non-perturbative behaviour. We can get a feeling for this behaviour 
if we make a simple ansatz for the hadronic spectral function, like e.g. the 
one predicted in the large-iY^ limit of QCD where the hadronic spectral 
function becomes an infinite sum of narrow states. [See e.g. Manohar's lectures 

in this school [2].] 

ilmn(i) = ^ f^M^5{t - Mf). (4.12) 
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It is easy to sec that at large values of a the function Tl{(j) goes to the mass 
squared of the ground state with the quantum numbers of the current J{x) 

nia) = ^'■'l \ Ml. (4.13) 

These considerations give us a feeling of what is the expected shape of 
the 7?.-function at asymptotic values of ct . It corresponds to the shape 
illustrated in Fig. 8. 




pQCD 



Nan~PaiturtMllvB 



Fig. 8. Expected shape of the short and long distance asymptotic behaviour of the 
function 7?.(<t). 



The problem of course is how to link the very short-distance regime (small 
cr) to the very long-distance (large a) one. We shall be able to say a 
little bit more about that after the discussion of non-perturbative power 
corrections in the next section. 

I wish to conclude this subsection with the discussion of an interesting 
model independent property of the 7?.-function [40]. It is the fact that, 
because of the positivity property of a spectral function ^lmli(t) > 0, 
the function — log A^(cr) is a concave function of cr; or in other words, the 
slope of the function 7?.(cr) must always be negative. This of course implies 
severe restrictions on the way that the two asymptotic regimes illustrated in 
Fig. 8 can be joined. The proof of this property is rather straightforward. 
It can be understood very simply by making an analogy with statistical 
mechanics: 'R-{a) can be viewed as the equilibrium "energy" {t) of a system 
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with variable "energy" t in thermal equilibrium with a second system at 
"temperature" 1/a. In this analogy, ^IvoHit) represents the "density of 
states" with "energy" t. Then the mean squared "energy fluctuation" is 
given by 

--^na) ^ {{t - {t)f) = it'') - {tf > 0, (4.14) 
which by definition is a positive quantity. 



4-2. Gaussian Transform Sum Rules 



The Gaussian transform of a spectral function 

g{s,T) = Hdie^p f-^^T^) -lnin(i), (4.15) 

V47rT Jo V 4t / TT 



can be evaluated in pQCD provided one keeps a finite width v2t resolu- 
tion in the gaussian kernel suSiciently large. These sum rules provide the 

framework to formulate quantitatively the notion of local duality. In the 
limit T where the gaussian kernel becomes a delta function 

^5{s-t), (4.16) 




the Gaussian transform coincides with the spectral function itself 

lim g{s, t) => -lmn(s). (4.17) 

There is in fact an interesting analogy between these sum rules and the 
theory of the "heat equation" [38]. The analogy is based on the observation 
that g{s,T) obeys the partial differential equation 

d'g{s,r) dg{s,r) 

which is the one dimensional heat equation if we interpret s as a pseudo- 
"position" variable and r as a pseudo-"time" variable. In this analogy, the 
hadronic spectral fimction corresponds to the initial "heat" (or "tempera- 
ture") distribution in a semi-infinite rod < s < oo and 5(s,t) measures 
the evolution in "time" r of the "heat" distribution in this rod. The cal- 
culation of g{s,T) in pQCD is an expansion in powers of Q!s(\/t) which 
blows up at small r-values; therefore, direct confrontation of the physical 
hadronic spectral function ;i:Imn(f) with pQCD cannot be made at small 
r-values. The confrontation can however be made at sufficiently large 
T-values by comparing the "evolved" physical hadronic spectral function 
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Q(s,t) with the corresponding pQCD theoretical prediction. The more we 
knew about non-perturbative corrections the smaller we should be able to 
go in "time" i.e. the smaller we would be able to take the width \/2t of 
the gaussian kernel and therefore the closer one could satisfy local duality. 

Let us discuss how to get the Gaussian transform from a generic two- 
point function like n(g^) in eq. (3.1). First, one evaluates n(g^) at a 
complex point = s + iA (s and A are real positive variables) and at its 
complex conjugate q-^ = s — iA (see Fig. 9) 







s+iA 

f 


Complex q 2 -Plane 


i 


< s-iA 



Fig. 9. Points in the complex g^-plane where the two-point function in eq. (4.19) is 
evaluated. 



and defines the combination, (I am assuming for simplicity that the disper- 
sion relation for Il{q'^) requires at most one subtraction, but the argument 
can be easily generalized as in the case discussed for the Laplace transform.) 

U(s + iA) U(s-iA) f°° , 1 1 

^ ' + ^ ' = / dt , -Im(t). (4.19) 



iA -iA Jo {t - S)2 -h A2 TT 

The integral in the r.h.s. brings in the convolution with a Lorentz-like 
kernel which we can write as a Laplace transform 



1 



{t - S)2 + A2 



CSO 



rfa;e-^^%-^(*-'*)'. (4.20) 



Applying to this integral representation the techniques developed in the 
previous section 4.1 allows us to construct the inverse Laplace transform 
operator which is needed to obtain the Gaussian transform in eq. (4.15) 
from the Lorentz transform in eq. (4.19). It is the operator 



lim lim 

iV— >(X) A^— >cx 



7T7^- (4-21) 



(AT-l)!^ ^ (5A2) 
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We then have the desired relation 



/47rr 



=2rL 



n(s + iA) ^ n(s - zA) 



iA 



/47rr 



dt exp I — 



-iA 



0^ 



4t 



(4.22) 



-lmn(i). 



I will later come back to the physics of Gaussian sum rules, but I need first 
to introduce some properties of finite energy sum rules which is the subject 
of the next section. 



4-3. Finite Energy Sum Rules 

The moments of spectral functions 

Mn{s) = [ ° dtf^ -lmn(i), (4.23) 
Ja 

with n = 0, 1, 2, • • • are also calculable quantities in QCD, provided that the 
upper limit so is sufficiently large. Finite energy sum rules follow from the 
analyticity properties of two-point functions already discussed in section 3. 
Applying Cauchy's integral formula to the contour in the complex g^-plane 
shown in Fig. 10 



Im q2 




Fig. 10. Contour of integration in the complex q'^-plane. 
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avoids the cut in the real axis and therefore results in 



(4.24) 



We now separate this integral into two pieces: one is the contribution from 
the paths above and below the real axis which pick up the discontinuity 
of the n(g^) function along this axis, and hence its imaginary part; the 
other is the contribution to the integral over the circle of radius sq [One 

should also check that there is no contribution from the little circle around the origin. 
The behaviour of n(g^) at small is usually known from low-energy properties of the 
two-point function.]. These two Contributions have to match, resulting in the 
equation 

/■So 1 r 



The integral in the l.h.s. is obtained using hadronic data, or inserting an 
hadronic ansatz; the integral in the r.h.s. is done in pQCD integrating 
the running coupling constant dependence over the circle of radius sq and 
scaling the result at /x^ = sq. Non-perturbative l/g^ -power corrections 
can also contribute in general to the r.h.s. integral. In QCD, and with d 
the asymptotic power behaviour of a given spectral function as defined in 
eq. (3.34), finite energy sum rules with n+d > 1 are particularly sensitive to 
non-perturbative power corrections. These corrections will be the subject 
of the next section. 

In finite energy sum rules, the question of subtractions in the primitive 
two-point functions can be dealt with by an appropriate integration by 
parts. For example, if one subtraction in n(g^) is required as for example 
in the case of the Adler function discussed in section 3.4, and we want to 
compute a general moment integral over the the circle of radius sq, we can 
use the identity {f{q^) is assumed to be non-singular,) 



with F{z) = Jq dwf{w). As one increases the n-power in a finite energy 
sum rule, one becomes more and more sensitive to the detailed high energy 
behaviour of the hadronic spectral function. In the QCD expression of 
the r.h.s., this is reflected by the appearence of higher and higher values 
of the coefficients in the perturbative series in as{s) and therefore worse 
convergence, as well as by the appearance of non-perturbative condensates 
of higher and higher dimension. 




(4.25) 




(4.26) 
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Finite energy sum rules provide a very useful way to do the correct 
matching between a given ansatz of the low-energy hadronic spectral func- 
tion and the onset of the perturbative continuum. For example, in the case 
of the Adler function discussed in section 3.4, let us consider the simplest 
hadronic ansatz of a narrow low-energy ground state plus a continuum i.e. 
the ansatz in eq. (2.23) which we have already discussed in connection with 
the phenomenological study of the Weinberg sum rules: 

llmni=i(i) = f^M'^5{t - Ml) + ^\0{t - so)[l + •••]. (4.27) 

In this case, the application of the finite energy sum rule in eq. (4.25) 
with n = and in the chiral limit where m„ = does not depend 

on non-perturbative condensates (because of the absence of gauge invari- 
ant operators of dimension d = 2 in QCD,) and results in the following 
constraint among the parameters of the hadronic ansatz 

/X'=2£^«o[l + •••]• (4.28) 

Numerically, using the experimental values of fp and Mp in this equation 

results in a value for the onset of the perturbative continuum sq ~ 1.8 GcV 
which is quite reasonable, in the sense that it is sufficiently large to trust 
pQCD [with gluonic corrections incorporated at the two loop level, the resulting so is 
reduced to sq ~ 1.6 GeV^.]. 

I can now comment on a very interesting connection between finite en- 
ergy sum rules and the Gaussian sum rules we discussed in the previous 
subsection. It has been shown in a general way [38] , that matching a given 
low-energy hadronic spectral function ansatz to the perturbative contin- 
uum via finite energy sum rules guarantees that the corresponding Gaus- 
sian transform at large widths has then the correct asymptotic behaviour 
predicted by pQCD. This is an important point in phenomenological ap- 
plications and therefore I want to discuss it in some detail. Very often 
one finds in the literature rather "precise" QCD sum rule "predictions" 
of masses or couplings which, however, have been obtained by fixing the 
onset of the pQCD continuum in the hadronic spectral function at some 
value using some "good physical sense" criteria, or even forgetting alto- 
gether about the QCD continuum. You can check that, unfortunately, the 
Gaussian transform of these spectral functions do not evolve to the shape 
predicted by pQCD. In other words, paramcterizations of spectral func- 
tions which are not correctly constrained by finite energy sum rules are not 
guaranteed to be globally dual to QCD and therefore should be avoided. 
Let us illustrate this with some pictures. Figure 11 shows the evolution 
in the pseudo-"time" variable r of the Gaussian transform of the spectral 
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function ansatz in eq. (4.27) with the onset of the continuum Sq fixed by 
the finite energy sum rule in eq. (4.28). In the "heat evolution" analogy 




Fig. 11. The Gaussian Transform of the Spectral Function in eq. (4.27) . 



the spectral function in eq. (4.27) corresponds to the initial "heat distri- 
bution" in the s-axis. The picture shows the evolution in "time" of this 
"heat distribution" in the interval 0.1 GeV^ < t < IGeV''. We observe 
that asymptotically in "time" , i.e. for r large, the spectral function evolves 
very well to the asymptotic "heat distribution" predicted by pQCD i.e. 

where erf(x) denotes the error function erf(a;) = dy . By con- 

trast, Fig. 12 shows the same evolution in the limit case of only a delta- 
function ansatz for the spectral function 

ilmni=i(i) = fpM^pS{t - Ml), (4.30) 

with no continuum. Clearly the corresponding asymptotic "heat distribu- 
tion" fails to reproduce the shape predicted by pQCD. Global duality of 
a given hadronic spectral function ansatz with QCD is only obtained pro- 
vided that the hadronic parameters are constrained to satisfy a system of 
finite energy sum rules equations. 
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Fig. 12. The Gaussian Transform of the Spectral Function in eq. (4.30) . 



5. Non Perturbative Power Corrections 

The physical vacuum of QCD is not the vacuum state which one uses 
in perturbation theory. Physical effects like spontaneous chiral symmetry 
breaking and/or confinement do not appear in an order by order pertur- 
bative treatment of QCD. A natural question to ask then is how are the 
pQCD results modified by non-perturbative effects at short-distances? This 
is the question which we shall be dealing with in this section. We shall see 
that non-perturbative effects in two-point functions evaluated at large Q^- 
values appear as inverse power corrections in Q^. They are the analogue 
of the higher twist corrections of deep inelastic scattering [2] . They can be 
systematically evaluated by doing Wilson's operator product expansion [41] 
(OPE) in the physical vacuum. The power corrections appear then as the 
product of Wilson coefficients, which are calculable perturbatively, times 
universal non-zero vacuum expectation values of gauge invariant operators, 
the so called condensates which, although excluded by definition order by 
order in perturbation theory, they can a priori have non-zero values in the 
physical vacuum. Typical examples are the lowest dimension quark con- 
densate i'tpip) and gluon condensate {tr Gf^^G^"). The operators which 
appear in any of the vacuum condensates which we shall be considering are 
assumed to be normal ordered. 
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5.1. Infrared Renormalon Ambiguities and Power Corrections 

The need of non-perturbative power corrections to colour singlet gauge 
invariant two-point functions can be "hinted" from what emerges already 
in perturbation theory when summing a specific class of contributions. Let 
us consider the calculation of a typical two-point function, e.g. the Adler 
function which we have already defined in section 3.4, and let us restrict 
our attention to the class of Feynman diagrams generated by the replace- 
ment of a virtual gluon by one effective charge chain as shown in Fig 13. 
Formally, the chain of bubbles in Fig. 13 corresponds to the replacement of 



Fig. 13. One of the Fcymiiaii diagrams which contributes to the Adler function with the 
gluon propagator replaced by an effective charge chain. 

the ordinary free gluon propagator by the full gluon propagator-like func- 
tion 



where ks denotes the euclidean virtual momentum carried by the chain; 
and ^ is the usual covariant gauge parameter. The interesting quantity 
here is the effective charge function aeff(^l;)- "^^^ existence of an effective 
charge in non-abelian gauge theories with properties precisely analogous to 
those of the well known effective charge of QED is discussed in refs. [42,43] 
and references therein. In terms of a characteristic scale, like e.g. the Ajg^ 
scale of the M S-renormalization scheme. 





(5.1) 




aeff(^l) = 



1 



(5.2) 




where 



/3i = 




(5.3) 
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The /?! factor is precisely the first coefiicicnt of the QCD /3-function. The 
effective charge encodes therefore the physics of the /3-function; in partic- 
ular the scale breaking property. Since, eventually, we are interested in the 
appearance of physical scales in QCD, it seems natural to focus our atten- 
tion on the properties of those diagrams of perturbation theory generated 
by the insertion of effective c/iar^e-exchanges. Of course, we do that at the 
simplest level of keeping only the one loop dependence of the /3-function. 
As far as one is only interested in general qualitative features this should 
not be a serious limitation. It is helpful now to look at the Q'^-k^ plot in 
Fig. 14, where denotes the euclidean energy momentum carried by the 
chain in Fig. 13. We shall refer to the integration regions < k% < Q"^ 




a1 

Fig. 14. Integration regions of the virtual euclidean momentum for a fixed euclidean 

external momenta Q^. 

and Q"^ < k'^ < oo shown in Fig. 14 as the IR region and the UV-rcgion 
respectively. It is this IR-region of integration which is at the origin of 
the appearance of singularities, the so called IR renormalons, when tree 
level gluon propagators are replaced by the full summation of perturbation 
theory diagrams which define a renormalon chain as illustrated in Fig. 13. 
For reference, we also show in the plot in Fig. 14 the lines = A| and 
— with Al = cxAj^ ~ 1 GeV the Landau scale at which the pQCD 
effective charge becomes singular. In pQCD calculations, the external eu- 
clidean momenta is always chosen to be Q'^ ^ A\; however, regardless of 
how large is taken, there will always be a region in the virtual k'^ inte- 
gration which is below A| where pQCD is not defined. It is precisely this 
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region of integration which, as we shaU see in the next paragraph, leads to 
IR renormalon poles in the so called Borel plane. The appearance of these 
singularities is welcome because they reflect the limitations of perturbation 
theory and indicate the presence of non-perturbative power corrections, as 
indeed the OPE in the physical vacuum suggests [The reader may be curious to 
know about the fate of the UV renormalons in QCD. Ref. [51] and the references therein 
discuss this interesting topic] . 

The leading contribution to the self-energy function ni=i(g^) in 
eq. (3.37) from the chain exchange diagrams in Fig. 13 in the IR domain of 
integration < fc|; < can be easily obtained and it has the following 

form (Cf = ^) 

N r'^^ 

ni=i(Q') ^ ^C^py^ dk%^a,s{k%). (5.4) 

Making the change of variables 

w/2 = -baa{Q^) log k%/Q^ , where bo = -/3i/27r (5.5) 
in the previous equation, one obtains 

n.-(e^)l,R = Jj-/°°.»e-^ ^ (^y. (5.6) 

which leads to a contribution to the Adler function 

AQ')\m = ^CF-^ Tdwe-^ -^(%) \ (5.7) 

This expression is already in the form of a Borel transform. An expansion 
around w = would generate the characteristic n! behaviour of the pertur- 
bative expansion in a(/i^) indicating that the corresponding series is not 
Borel summable. As already discussed by other authors [44-46], we find 
that the leading IR renormalon contribution to the Adler function appears 
as a pole in the Borel plane at the location w = 2. There is no term in the 
IR expansion of the domain of integration which leads to a pole at w = 1. 
From the previous change of variables one also sees that low values of w are 
associated with momenta around the large scale Q^, where perturbation 
theory is expected to give a good description of the dynamics. However, 
as w goes up (and it has to go all the way up to infinity) one enters deeper 
and deeper into the IR region, where perturbation theory must fail or else 
QCD would not describe the spectrum of hadronic bound states that are 
observed. The singularity at w = 2 exhibits this in its crudest form. More 
precisely, this singularity implies an ambiguity in the perturbative evalu- 
ation of the Adler function. Therefore, the analytic continuation in w in 
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eq. (5.7), or cquivalently the resummation of perturbation theory into the 
effective charge agg(fc|;) of eq. (5.2), that one has tried in order to obtain 
the full solution has failed. The ambiguity is encoded in the unavoidable 
prescription to skip the pole; and is of the form <^ Q-'^/boaiQ) . 

An(g)|iR^^-^^^-^(^-J , (5.8) 

i.e., an ambiguity which has the same ^ pattern as the ghion conden- 
sate contribution which appears in the OPE evaluation of the Adler func- 
tion [47-50]. Since the Adler function must be an unambiguous physical 
observable, there must exist another contribution that cancels the one in 
eq. (5.8). In other words, (Borel) resummed perturbation theory is requir- 
ing the presence of the gluon condensate (which is also ambiguous for the 
same reason) to combine with the result of eq. (5.7) and yield a final well- 
defined answer [44]. This is an example of how all-orders perturbation 
theory is capable of "hinting" at nonperturbative dynamics. 



5.2. The ITEP Power Corrections 



The intuitive physical picture which lies behind the non-perturbative power 
corrections introduced by the ITEP group [37] can be understood as follows. 
In terms of Feynman diagrams, the configurations which are sensitive to the 
IR-region are those where the hard external momentum Q flows into the 
diagram and goes along the quark loop with very little transfer to the gluon 
propagator, or the one-renormalon chain. The gluon propagator is then in 
a non-perturbative regime where many soft gluon interactions are possible 
since the coupling constant in this regime is very large. The idea of the OPE 
is to factor out the hard part of the Q-flow and calculate it perturbatively, 
while collecting the bulk of the soft non pcrturbativc gluon interactions in 
a condensate of gluons which will be treated as a phenomenological param- 
eter. In practice, the calculations in the OPE are performed by integrating 
the fermion loop in the presence of "background" gluon fields, including 
hard perturbative gluonic corrections if necessary [Ref.[52] gives a lot of tech- 
nical details on how to do these calculations.] . The average of all possible gauge 
invariant configurations of "background" gluon fields defines then the con- 
densates. Background fields in an ordinary Feynman diagram are indicated 
by terminating dots, like in Fig. 15. The dots are a short-hand notation 
for the non-perturbative average over soft interactions. The perturbative 
evaluation of the "hard" part of the loop gives the corresponding Wilson 

coefficient and the net result is a contribution of the type Cgg— — q4 , 
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Fig. 15. Feynman diagram representation of gluon condensate contributions. 

where Cqg denotes the Wilson coefficient. The Wilson coefficient depends 
on the specific two-point function one is calculating, while the gluon con- 
densate g'tiGa'^ Gl^^j) is treated as a universal phenomenological parameter. 
For example, in the case of the Adler function discussed in section 3.4 the 
gluon condensate gives a correction to the perturbative result in eq. (3.40): 



A{Q)\aG = Y^- . (5.9) 



For two-point functions with composite operators of light quarks, there 
are also l/Q^-corrections induced by the quark condensate, as illustrated 
in Fig. 16 below which give contributions of the type C^^ , 




Fig. 16. Feynman diagram representation of quark condensate contributions. 

where, as before, the Wilson coefficient C^^ depends on the specific current 
operator. The appearance of a mass factor in front of the quark condensate 
is due to the conservation of chirality. In the case of the Adler function 
discussed in section 3.4 the quark condensate gives a correction to the per- 
turbative result in eq. (3.40): 

mujuu) + md{ 

Higher order terms in the OPE bring in contributions from condensates 

of higher and higher dimension. Three new types of condensates appear 
at 0(1/0^): the quark gl uon mixed condensate illustrated in Fig. 17 be- 
low which brings contributions of the type — — ^r^jj^.^ ^j^g 



AQ')U^ = ' ' ■ (5.10) 
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Fig. 17. Feynman diagram representation of contributions from the mixed quark gluon 
condensate. 



three gluon condensate illustrated in Fig. 18 which can give contributions 




Fig. 18. Feynman diagram representation of contributions from the three gluon conden- 
sate. 



like Cqgg — ^ — QB ; and the four quark condensate illustrated in 

Fig. 19 which contributes to terms like C^^^^ {^r^n^rip) ^ ^j^^ 




+ 



Fig. 19. Feynman diagram representation of contributions from the four quark conden- 



Adler function discussed in section 3.4 the quark gluon mixed condensate 

2 

is suppressed by ^ factors. It has been shown [53,54] that the three gluon 
condensate, without extra hard gluon corrections, gives no contribution to 
the Adler function. The four quark condensate contribution was first cal- 
culated in ref.[37] and it gives a correction to the perturbative result in 
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eq. (3.40): 



3 / 

2 \(w7/:i75A''u - rf7^75A"rf)' 



\ q=u,d,s I 



(5.11) 



To leading order in the l/N^ expansion these four quark condensates fac- 
torize into products of the lowest dimension quark condensate. 

The phenomenological determinations of the gluon condensate and the 
quark condensate have still rather big errors. A generous range for the 
gluon condensate which covers most of the determinations is 

(a,GrG;iJ = (0.08 ± 0.04) GcV^ (5.12) 

The quark condensate is scale dependent. The values obtained from a 
variety of sum rules are within the range 

{uu + dd){lGeV') ^ _(o.oi3± 0.004) GeV^ (5.13) 

It turns out that with those values for the lower dimension condensates, 
the contribution from the leading non-perturbative power corrections in 

most of the sum rule applications -at the values where one can trust 
the pQCD contributions- turn out to be already relatively small. This is 
due to the fact that the scale of perturbation theory is rather large; 
therefore, in order to get a convergent series in as(Q^) one is obliged to go 
to rather high values of where the power corrections are already quite 
small. 

5.3. The Weinberg Sum Rules in the Light of the OPE 

We are now in the position to understand why it is that the Weinberg sum 
rules we discussed in section 2.3 are satisfied in QCD. For this it is sufficient 
to consider the OPE applied to the two-point function YlLR{q) defined in 

eq. (2.12) [Notice that the OPE of this self— energy function is free from rcnormalon 
ambiguities, because in perturbation theory it is protected by chiral symmetry, which 
ensures that niji(Q^) = order by order in powers of a^.]. Since this two— point 
function is chirally symmetric when quark masses are set to zero, the only 
possible inverse powers of which can appear are those modulated by lo- 
cal order parameters of S^SB. In QCD there are no local order parameters 
of dimensions d = 2 and d = 4. Notice that (^/'V'), which has dimensions 
d = 3 contributes to ^-terms but it appears multiplied by a quark mass 
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factor; and therefore it disappears in the chiral limit. The first possible 
contribution conies from local order parameters of dimension d = 6, i.e. 
the four-quark condensate contributions of the type illustrated in Fig. 18. 
The specific contribution to H-LniQ^) has been calculated in ref. [37]. Their 
result is particularly simple in the large- A^c limit where the vacuum expec- 
tation value of four-quark operators factorizes in a product of (V'V') terms, 
and therefore we have that 

lim^n.,(Q^) = ' [-4.^ + OiaD) i^^f] ^ O (') . (5.14) 



From this follow the two constraints in eqs. (2.18) and (2.19), and hence 
the two Weinberg sum rules in eqs. (2.15) and (2.16). 

We also observe that the YIlr{Q^) expression in eq. (2.33) which re- 
sults from the resonance dominance ansatz for the vector and axial-vector 
spectral functions has the required behaviour at large-Q^ 

^lirn^ n^H(Q^) = (/^M« - flM\) l^ + 0{^). (5.15) 

This is not surprising since the ansatz in question is nothing but a specific 
choice compatible with the spectrum predicted by QCD in the large-A^c 
limit. 

In section 2.3 I stressed the fact that it was not enough to know the 
the very low behaviour of li-LRiQ"^) in order to calculate the integral in 
eq. (2.32) which gives the electromagnetic 7r+ — 7r° mass difFerence. We now 
observe that it is not enough cither to know only the very high behaviour 
of the same function. The ^ behaviour of the short-distance OPE in 
eq. (5.14) cannot be extrapolated to the infrared region of the integral 
since it would lead to a divergent result as well. The simple resonance 
dominance ansatz in eq. (2.33) provides however, in this case, a correct 
matching of short-distance behaviour with the long-distance behaviour. 



6. Some Examples of QCD Sum Rules 

Here we shall discuss two applications of the techniques about sum rules in 

QCD developed in the previous sections. From a phcnomenological point 
of view, the most striking application is the determination of the QCD 
coupling constant as(m^) which has been made using the LEP data on 
hadronic tau decays. Another recent application uses very general proper- 
ties of QCD sum rules to set lower bounds on the light quark masses. 
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6.1. Determination of the QCD Coupling Constant from Tau decays 

Hadronic tau decays are governed by the spectral function associated with 

the charged current L^{x) = u{x)^^^^{\ — 75)de(x), where dg denotes the 
Cabibbo rotated d-quark. From the QCD point of view, the relevant two- 
point function is then 

nril?) =4i / rf^a;e'«-(0|T(L^(x)L'^(0)t) |0) 

J (6.1) 

which in full generality, when the quark masses are not neglected has two 
invariant functions corresponding to physical states with total angular mo- 
mentum J — \ and J = 0. From the phcnomcnological point of view, the 
quantity which is accessible to experimental determination is the hadronic 
tau decay branching ratio 



_ T{T^Vr+ hadrons) 

r(T Vre Ve) 







(6.2) 



The basic idea to apply here is the one we explained in the section 4.3 on 
finite energy sum rules, i.e. to confront the experimental determination of 
Rr with the corresponding theoretical result [55] obtained from the contour 
integral in the complex g^-plane shown in Fig. 10 with sq = rn^: 



2 ^ 2 



(6.3) 



The analyticity properties of the two-point function in eq. (6.1) imply 
that the r.h.s.'s of cqs. (6.2) and (6.3) should be equal. Notice that the 
evaluation of this integral only requires the knowledge of the invariant func- 
tions n^^^((7^) and I\!"j^\[q^) in the complex big circle of radius \q^\ = rn^ 
and vrir — 1.78 GeV is a reasonably large scale to apply pQCD. Non- 
perturbative power corrections to this integral bring inverse powers of m^. 
In fact, one of the nice features about this integral is that because of the 
combination of g^-powers in the phase space factor which modulate the 
ni,z,(g^)'s functions in the integrand, the leading 1/m^ non-perturbative 



power corrections are suppressed up to terms of O 



^ {asGG) 
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Also the contribution from the region in the circle near the real axis, where 
the validity of the OPE is perhaps questionable, is suppressed by the kine- 
matic factor ^1 — . The result of the explicit calculation, which in- 
vokes quite a lot of refined techniques [56] , can be written as follows 

Rr = iVc (IKdP + \Vus?f ^EW {l + <5ew + '5^°^ + <5np} , (6.4) 



where 5'ew = 1.0194 and 5-^-^ = 0.0010 are the contributions from the 
leading and next-to-leading electroweak corrections, and 

=^^lS!!lri + 5.2023 ( i^^^ ) ' + 26.366 



O 



as(m^) 



(6.5) 



is the result of the pQCD calculation in the chiral limit. The remaining 
factor S-^p w —0.016 ± 0.005 includes the estimated effects of small quark 
mass-corrections and non perturbative power corrections. 

The most recent analysis, using the ALEPH data [57] gives the value 

as{nir) = 0.348 ± O.OOScxp ± O.OlSui- (6.6) 

This value, when evolved with the renormalization group running to the 
Z-mass scale gives 

as{Mz) = 0.1211 ± O.OOOSexp ± 0.0016th ± O.OOlOevoi, (6.7) 

which, with these errors, represents the most accurate determination of 
as{Mz) from a single experiment . 

6.2. Lower Bounds on the Light Quark Masses 

The light u, d and s quarks in the QCD Lagrangian have bilinear couplings 
-Cqcd = C^Qcn - X! '^q^i^)lix) , (6.8) 

q=u,d,s 

with masses ^ which explicitly break the chiral symmetry of the La- 
grangian. The present situation concerning the values of the light quark 
masses in the MS renormalization scheme, or any other scheme for that 
matter, has become rather confusing because of recent determinations of 
the light quark masses from lattice QCD [58,59] which find substantially 
lower values than those previously obtained using a variety of QCD sum 
rules, (sec e.g. refs. [60-62] and earlier references therein.) These new lat- 
tice determinations are also in disagreement with earlier lattice results [63] . 
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More recently, an independent lattice QCD determination [64] using dy- 
namical Wilson fermions finds results which, for m„ + are in agreement 
within errors with those of refs. [58,59]; while for they agree rather 
well, again within errors, with the QCD sum rule determinations. [For a 
very lucid introduction to the techniques of Lattice QCD, see the lectures of M. Liischer 
at this school [65].] Here, and as an application of the topics discussed in the 
previous sections, I shall show that there exist rigorous lower bounds on 
how small the light quark masses rus + niu and + mu can be [66] . There 
are two-point functions 

= i J d'^xe'^-^0\T {0{x)0{Q)^) |0), (6.9) 

with an appropriate choice of the local operator 0(x) which are particu- 
larly sensitive to quark masses: the case where O is the divergence of the 
strangeness changing axial current 

0{x) = ^^A^'{x) = (m, + m„) :s{x)i-f5u{x):; (6.10) 

and another one where O is the scalar operator S{x), defined as the isos- 
inglet component of the mass term in the QCD Lagrangian 

0{x) = S{x) = m\:u{x)u{x): + :d{x)d{x):], (6.11) 

with 

m = ^(m„ + md). (6-12) 

Both two-point functions obey dispersion relations which in QCD require 
two subtractions, and it is therefore appropriate to consider their second 
derivatives (Q^ = —q^)- 

2 1 

The bounds follow from the restriction of the sum over all possible hadronic 
states which can contribute to the spectral function to the state(s) with the 
lowest invariant mass. It turns out that, for the two operators in (6.10) 
and (6.11), these hadronic contributions are well known phenomenologi- 
cally; either from experiment or from chiral perturbation theory (xPT) 
calculations. On the QCD side of the dispersion relation, the two-point 
functions in question where the quark masses appear as an overall fac- 
tor, are known in the deep euclidean region from pQCD at the four loop 
level. The leading non-perturbative ^-power corrections which appear 
in the operator product expansion when the time ordered product in (6.9) 
is evaluated in the physical vacuum [37] are also known. For the sake of 



= / dt ^^ , -lm^(t). (6.13) 
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simplicity, I shall only discuss here the case where O is the divergence of 
the strangeness changing axial current in eq. (6.10). 

We shall call ^'5(0'^) the two-point function in (6.9) associated with the 
divergence of the strangeness changing axial current in (6.10). The lowest 
hadronic state which contributes to the corresponding spectral function is 
the ii'-pole. From eq. (6.13) we then have 



2 



(Mi + 02)3 



VkM^ + / 

Jto 



dt 



1 



-Im*5(i), (6.14) 



where to — {Mk + 2171^,-)'^ is the threshold of the hadronic continuum. 

It is convenient to introduce the moments 5]jv((5^) of the hadronic con- 
tinuum integral 



Siv(O') 



f 

Jto 



dt 



to + Q^' 

{t + 02)3 \ t + Q^ 



N 



-Im*5(t). 

TT 



(6.15) 



One is then confronted with a typical moment problem (sec e.g. rcf. [67].) 
The positivity of the continuum spectral function ilm5'5(t) constrains the 
moments Ejv(Q2) and hence the l.h.s. of (6.14) where the light quark 
masses appear. The most general constraints among the first three mo- 
ments for = 0, 1, 2 are: 



So(O^) > 0, Si(O') > 0, S2(0^) > 0; 



(6.16) 



Eo(O') - Si(O') > 0, Si(O') - S2(0') > 0; 



(6.17) 



So(0')S2(0')-(Si(o'))'>o. 



(6.18) 



The inequalities in eq. (6.17) are in fact trivial unless 20^ < ^o, which 
constrains the region in to too small values for pQCD to be applicable. 
The other inequalities lead however to interesting bounds which we next 
discuss. 

The inequality T,o{Q^) > results in a first bound on the running 
masses: 

[m.(0^)+m„(0^)]^>^^ 

(6.19) 



1 



1 



1 



Ml 



1 I 11 Q»(Q^) , 

^ 3 TT 



where the dots represent higher order terms which have been calculated 
up to 0{ag), as well as non-perturbative power corrections of O and 
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strange quark mass corrections of O and O (^^^ including 0{as) 

terms . Notice that this bound is non-trivial in the large -N^ limit (/|. ^ 
0{Nc)) and in the chiral limit (m^ ~ ^k)- The bound is of course a 
function of the choice of the euclidean Q-value at which the r.h.s. in 
eq. (6.19) is evaluated. For the bound to be meaningful, the choice of Q 
has to be made sufficiently large. In ref. [66] it is shown that Q > 1.4 GeV 
is already a safe choice to trust the pQCD corrections as such. The lower 
bound which follows from eq. (6.19) for m,j + m., at a renormalization scale 
= 4 GeV'^ results in the solid curves shown in Fig. 20 below. These 



120 
> 100 




1.4 1.6 1.8 2 2.2 2.4 

Q[GeV] 

Fig. 20. Lower bound in MeV for [ms +m„](4GeV^) versus Q(GeV) from eq. (6.19). 



are the lower bounds obtained by letting vary [68] between 290 MeV 
(the upper curve) and 380 Mev (the lower curve). Values of the quark 
masses below the solid curves in Fig. 20 are forbidden by the bounds. For 
comparison, the horizontal lines in the figure correspond to the central 
values obtained by the authors of ref. [58]: their "quenched" result is the 
horizontal upper line; their "unquenched" result the horizontal lower line. 
It must be emphasized that in Fig. 20 one is comparing what is meant to 
be a "calculation" of the quark masses -the horizontal lines which are the 
lattice results of ref. [58]- with a bound which in fact can only be saturated 
in the limit where So((5^) = 0. Physically, this limit corresponds to the 
extreme case where the hadronic spectral function from the continuum 
of states is totally neglected! What the plots show is that, even in that 
extreme limit, and for values of Q in the range 1.4 GeV < Q < 1.6 GeV the 
"unquenched" lattice results of ref. [58] are already in difficulty with the 
bounds. 

The quadratic inequality in (6.18) results in improved lower boimds for 
the quark masses which we show in Fig. 21 below. Like in Fig. 20, these 
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are the lower bounds obtained by letting vary [68] between 290 MeV 

(the upper curve) and 380 Mev (the lower curve). Values of the quark 
masses below the solid curves in Fig. 21 are forbidden by the bounds. The 
horizontal lines in this figure are also the same lattice results as in Fig. 20. 

The quadratic bound is saturated for a (5-like spectral function repre- 
sentation of the hadronic continuum of states at an arbitrary position and 
with an arbitrary weight. This is certainly less restrictive than the extreme 
limit with the full hadronic continuum neglected, and it is therefore not 
surprising that the quadratic bound happens to be better than the ones 
from T,m{Q'^) for iV = 0, 1, and 2. Notice however that the quadratic bound 
in Fig. 21 is plotted at higher Q-values than the bound in Fig. 20. This 
is due to the fact that the coefficients of the perturbative series in asiQ"^) 
become larger for the higher moments. In rcf [66] it is shown that for the 
evaluation of the quadratic bound Q > 2 GeV is already a safe choice. One 
finds that in this case even the quenched lattice results of refs. [58,59] are 
in difficulty with these bounds. 

Similar bounds can be obtained for rriu + rud when one considers the 
two-point function associated with the divergence of the axial current 

d^A^ix) = (md + TO„) :J(x)«75w(,t): . (6.20) 

The method to derive the bounds is exactly the same as the one discussed 
above and therefore we only show, in Fig. 22 below, the results for the cor- 
responding lower bounds which one obtains from the quadratic inequality. 
One finds again that the lattice QCD results of refs. [58,59,64] for niu + md 
are in serious difiiculties with these bounds. The lower bounds discussed 
above are perfectly compatible with the sum rules results of refs. [60-62] 
and earlier references therein. 
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1.4 1.6 1.8 2 2.2 2.4 

Q[GeV] 

Fig. 22. Lower bound in MeV for [md + mti](4 GeV"^) from the quadratic inequality. 

The two examples discussed in this section are rather illustrative of 
how sum rules can relate QCD short distance calculations to low energy 
hadronic experimental information. There are many other examples dis- 
cussed in the literature. My hope is that these introductory lectures may 
help you to read the literature and to find new applications as well. 
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